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' Flux-pinning-induced stress and strain distributions in a thin disk superconductor in a perpen- 

I dicular magnetic field is analyzed. We calculate the body forces, solve the magneto-elastic problem 

■ and derive formulas for all stress and strain components, including the magnetostriction A_R/_R. The 
flux and current density profiles in the disk are assumed to follow the Bean model. During a cycle 
of the applied field the maximum tensile stress is found to occur approximately midway between 
the maximum field and the remanent state. An effective relationship between this overall maximum 
stress and the peak field is found. 

pH : I. INTRODUCTION 

o 

^ During recent years top-seeded growth of bulk Y-Ba2Cu30i/ and (RE)-Ba2Cu30i/ (RE — rare earth elements) 

5h has developed enormously towards production of large single-grain superconductors with strong flux pinning^^^. An 
impressive critical current density value of Jc — 100 kA/cm'^ at 77 K has now been obtained for bulk (Nd-Sni-Gd)- 
^ Ba2Cu30y with 10 mol% of Gd-211.^ Large-grain superconductors of this type can trap very high magnetic fields 
I t _ in the remanent state, exceeding the field from conventional permanent magnets by an order of magnitude and even 
I more. Therefore, it is today strongly believed that the field-trapping ability of bulk superconductors is useful in many 
^ ' practical applications like magnetic separators, magnetron sputtering systems, motors etc. 

A major problem of such bulks is that mechanical fracturing frequently occurs during the magnetic activation 
process. The strong flux pinning unavoidably leads to large body forces, which during descent of the applied field are 
tensile. Due to the irreversible nature of the flux dynamics, the resulting stress distribution becomes non-uniform, 
Q ■ and can locally reach the tensile strength of the material so that fatal cracking takes place. A systematic study of 
this problem was first made by Ren et al.'^, who also made model calculations of the internal stress in a cylindrical 
superconductor during the field descent after field-cooling. Later the modelling of the magneto-elastic behaviour has 
been extended considerably by including various geometrical cases and magnetization conditions^^^. 
. Common to all previous modelling work is that the superconductor is approximated as infinitely long in the direction 
' of the applied field. Today, one sees that single-grain disks with diameters up to 6 cm are grown on routine basis, 
whereas low- gravity environment has allowed growth of even larger disks - up to 12.7 cm in diameter and 2 cm in 
C thickness. When the diameter becomes much larger than the thickness the infinitely long cylinder approximation 

^ is known to give a poor description of fiux and current distributions. A better approximation can then be a very thin 

■ sample in perpendicular field, where exact analytical results for flux distribution are also available. They were used 
' to calculate magnetostriction for thin long strips in Ref. 14, but only for the virgin branch. In this work we use the 

■ thin disk solution, and flnd all the stresses and strains for the full magnetization cycle. We present here results of this 
' model and discuss the relationship between the peak applied field and the overall maximum tensile stress occurring 

^ in the disk during pulsed field magnetic activation. 

O . II. MAGNETO-ELASTIC PROBLEM 
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Consider a thin superconducting circular disk of radius R and thickness d, where d R. The superconductor is 
placed in a magnetic field, Ba, applied perpendicular to the disk plane, i.e., along the z-axis. When the local field at the 
5^ • edge exceeds the lower critical field, B d, the vortices start to penetrate into the disk. Due to strong demagnetization 
d I effects this happens already at Ba > \/ djR Bci ■ The presence of pinning centers will cause the vortices to distribute 
non-uniformly over the disk volume. According to the critical state model the quasi-stationary flux distribution 
adjusts irreversibly so that the local current density never exceeds the critical magnitude Jc- In analyses of this 
behaviour a common simpliflcation is to ignore the effects of a lower critical field, and also the presence of any 
reversible magnetization. Under these conditions analytical results for the field and current distributions have been 
derived for the thin disk geometry. 
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FIG. 1. Internal stresses and body forces in a superconducting disk placed in a perpendicular applied magnetic field. 



The force per unit volume acting on the vortices is given by the Lorentz force f = J x B, where J and B are the 
local current and flux density, respectively. In quasi-static equilibrium this force is balanced by pinning forces, and 
hence, f also represents the distribution of body forces transmitted onto the crystal lattice in a given magnetized 
state. The symmetry of the disk geometry implies that the induced current has only an azimuthal component, and 
the force averaged over the disk thickness points in the radial direction, see Fig. 1. 

For the thin disk one can therefore write 



/(r) - J{r)B{r), (1) 

where J{r) is the current density averaged over d, and B{r) = Bz{z — 0) is the flux density in the disk's center plane. 

Exact analytical results for a thin disk in perpendicular field are available only for the Bean model, where Jc is 
assumed field- independent. We adopt this model in the present analysis, and review here only the results needed for 
the calculations of the magneto-elastic behavior. 

Starting with a fiux free superconductor the application of a magnetic field, Ba^ will cause penetration of flux into 
the disk up to the point r — a, where 

R 

and Bd — (2) 



COSh{Ba/Bd 

In the annular vortex-filled region the fiux density has the strongly non-linear profile given by 

B(r) Ba w., , ^ 2 l-Ocote 

^ ' - arccosh(R/r) -I- - / — — dO, a < r < R. (3) 



Bd Bd TT Jarcsin(a/r) a/ 1 - (r/R)^ sin 6* 

The central region r < a has -B = 0, as for the long cylinder. 

The shielding current is for a thin disk distributed over the entire area according to 




Jir)/Jc=l 2 (r lR^^a^\ (4) 



Note here that J is not proportional to the flux density gradient dB/dr, thus preventing us from using the relation 
fdr = — c?(i3^)/2/io, which is valid for a long cylinder. Figure 2 (top) shows graphs of the body force distribution, 
Eq.(l), during virgin magnetization of the disk. Since the current has a constant magnitude where the flux density is 
non-zero, the force profile has the same shape as the function B{r). 

When the applied field is reduced after reaching a maximum value B'^^^ , the disk undergoes remagnetization where 
the front is located at r = 6, with 

" cosh - Ba)/2Bd] ■ ^ ' 

In the inner region, r < 5, the fiux density generated by the maximum field is unchanged, whereas in the outer 
annulus, b < r < R, considerable redistribution of fiux is taking place. The current, on the other hand, becomes 
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modified in the entire disk area. During field reduction the disk is tlierefore divided into three regions, each described 
by its own set of functions J(r), B{r) and thereby also /(r). One has for: 
Region I (outer); b < r < R; 



:^ = §i+arccosh(R/r) + - 



Region II (middle); a <r < b; 

B{r) B^^'^ 



arcsin(b/r) 



1 - 61 cot 6* 



arcsmf a, 



/r) ^Jl-{v/Rfsme 



-I 

o ai 



r/2 



1 -Qzoid 



■Ad 



i(b/r) a/1 - (r/R)2sin 

J(r)/Jc = 1 



(6) 



Region III (inner); r < a; 



arccosh(R/r) + — / 

J(r) 




4 / r 

= — 1 H — arctan | —\ 

TT 



(7) 



B{r) = 0, 



J(r) 

Jc 



2 arctan — 




arctan — 




(8) 



Shown in Fig. 2 (bottom) arc plots of /(r) at various stages of field reduction. A cusped peak in the body force 
is seen to accompany the remagnctization front at r = 6. After reversing the field sweep the outer part of the disk 
immediately starts to experience forces that will create tension. When the field has been reduced to zero considerable 
remanent body forces, a mixture of tensile and compressive forces, remain in the disk. 

To determine the quantitative magnetostrictive behaviour a separate analysis of the thin disk elasticity problem is 
required. For simplicity, we assume that the superconductor is mechanically isotropic, and that the strains are well 
below the fracture limit allowing linear elasticity theory to be applicable. As mechanical boundary conditions for the 
disk we choose a free surface, i.e., the normal stresses, ar and az vanish at the surfaces. Moreover, since the disk 
is thin one can assume that cr^ = throughout the volume, and that the deformation is described by only a radial 
displacement field, u{r). The strain components are then 



Br = u'{r) and eg = u/r, 
which are related to the non- vanishing stresses cFj. and (Tq by 

Ecr = Or — voo and Ee^ = ao — vcTt 



(9) 



(10) 



Here E and v is the Young's modulus and Poisson's ratio, respectively. The condition of static equilibrium, see Fig. 1, 
is that 



CTr — ere 



+ / = 0, 



which in terms of the displacement field is expressed as 



u" + -u' u- 



E 



-/ = o. 



or 



d_ 
dr 



1 d{ur) 
r dr 



E 



(11) 
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(13) 
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FIG. 2. Body force distribution in a disk during field rise from the virgin state (top), and subsequent field descent (bottom). 
Positive force means that it points radially out from the center. 



Integrating twice, and using that arir = R) = 0, and u(0) — 0, one obtains 



1-1/ l + u 



'2 J- J ' 

2 ; / r fdr 



It immediately follows that the external dilatation Ai? ~ u{R) is given by 

AR l-u I 



R E R^ Jo 



'fdr. 



The internal strains can be expressed as 



er{r) _ l-v 
ee{r) ^ 2E 



(14) 



(15) 



(16) 



with the upper and lower sign corresponding to and eg, respectively. Similarly, the radial stress and the hoop stress 
are given by 



(Tr{r) _ 1 ~ V 

o-e(r) ~ 2 



(17) 
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The eqs. (14-17) arc exact expressions from which the magneto-clastic response can be calculated once the force 
distribution /(r) is known. One can see immediately that the two strains and the two stresses are equal in the center, 
i.e., er(0) = ee{0) and (7^(0) = cre(O), independent of the magnetized state. Notice also that the hoop stress at the rim 
of the disk is proportional to the dilatation, ae^R) = EAR/R, as it should according to Eq. (10) and the condition 
ar{R) = 0. 



III. RESULTS AND DISCUSSION 



A. Magnetostriction AR/R 

The irreversible behaviour of the magnetostriction, AR/R, as a function of the applied field, Ba, is calculated from 
Eq. (15) using the expressions for J(r) and B{r) that were listed in the previous section. The result is shown in Fig. 3 
(left), where hysteresis loops with three different maximum fields are plotted. The overall shape of the loops is quite 
similar to that found in the long cylinder casc.^ When the maximum applied field becomes much higher than Ed a 
large part of the loop is well represented by two straight lines forming a cross centered a.t Ba = 0. Asymptotically, 
this linear behaviour is described by 

^ = l^..«S.(.^ + 7.57.10-»). (18, 

Therefore, the width of large field {B^'^^ » Bd) hysteresis loop rapidly approaches 

{AR/R)^ - {AR/R)^ = 2^ J,RB,. (19) 

Interestingly, the same simple relation holds also for the long cylinder when S™*^ > 2Bp, where Bp = fxoJcR is the 
full penetration field. ^ The virgin branches of the loops are not included in Fig. 3 since on the chosen scale they are 
barely resolved. 

The corresponding result for thin long strips is given by Eqs. 16-17 of Ref. 14. For large B™^^ it reduces to our 
Eq. (19), except that the factor 1 — v is missing. This factor should appear when the stress from the currents on the 
far edges of the strip is taken into account, as also discussed in Ref. 14. 

Dilatation loops found experimentally-'^^'"'^^ can deviate from the shown Fig. 3 (left) due to a field dependence of 
the critical current. If this dependence is rather strong, the loops will narrow down for large fields, as shown in 
Fig. 3 (right). These curves were calculated assuming an exponential dependence of the critical current, Jc{B) = 
Jc(0) exp(— |i?|/2i?d)- Flux and current distributions cannot be written down explicitly in case of field-dependent Jc, 
however they can be found numerically by solving a set of integral equations given in Ref. 15. It allows calculation 
of magnetization loops, as explained in Ref. 16, where M decreases with B at large fields in accordance with chosen 
Jc{B) dependence. A similar procedure was used to calculate the dilatation loops of Fig. 3 (right). 
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B. Stress Distributions 



As the applied field increases, the body force points towards the disk centre, leading to stresses that everywhere 
are compressive (negative). During the ramp up of the field the compressive stress grows monotonically, and always 
with a spatial distribution where the central unpenetrated region experiences the largest compression. Since these 
stress values are normally much lower than the compressive strength of the superconductor, we focus here instead on 
the descending field stage of the magnetic activation process. 




As an example, let us consider the case where the maximum applied field has the value B^^'^^ = ABd- When Ba 
starts to decrease, the stress distribution becomes more complex since the body force turns expansive in the outer 
remagnetized region. The Fig. 4 illustrates how the stresses develop during the field descent. Very soon a peak of 
tensile (positive) stress appears near the remagntization front, and the peak value grows in height. At an intermediate 
field the peak tension reaches an overall maximum, which for the radial stress amounts to 0.8ao, where ctq = BdJcR- 
The maximum hoop stress is slightly smaller. It is during this part of the field descent the disk is most likely to crack, 
namely the first time the material encounters a stress equal to the tensile strength. 
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FIG. 5. Time evolution of maximal tensile stress in the 
disk during a full cycle of applied field starting from the vir- 
gin state. 



FIG. 6. The overall maximal stress reached during the full 
field cycle, as found from Fig. 5, versus the applied field am- 
plitude. The dimensional values on the right and top axes are 
given for a typical disk with parameters given in the text. 



Shown in Fig. 5 is the time evolution of the stresses during one full cycle of the applied field. The curves show 
the maximum tension occurring in the disk. The parts given zero value correspond to the stages where the stress is 
everywhere compressive. It becomes evident that the most 'dangerous' stage of the process is slightly after midway 
between the maximum field and the remanent state. The tension in the remanent state is roughly twice smaller than 
the overall maximum value, which we denote as cr^^P^'^'^i^^ _ The next Figure shows how this cr^"?'^''™'*'^ depends on the 
maximum applied field in the cycle. The dependence is essentially parabolic, and is well fitted by a simple formula. 



BdJcR 



0.063 



O.OSE 



(20) 



By substituting the following realistic numbers for disk parameters; = 10 kA/cm^, d = 2 cm, i? = 3 cm, one 
obtains dimensional values for the peak tension versus the maximal applied field as shown in the top and right axes of 
Fig. 6. This graph represents the essence of the present work, and may hopefully serve as a useful guide in practical 
situations. E.g. it follows that for a disk with tensile strength of 20 MPa,^^ one has to apply a field less than 18 T to 
avoid cracking during the pulsed-field activation. 

The results presented here are obtained in the limit of thin disk in perpendicular field. Combined with the previous 
results for a long cyliner^, they give two extreme cases between which all real samples fall into. Unfortunately, 
the analytical treatment for samples with finite thickness is not possible, and numerical results are not available. 
Therefore, we may only hope that the key quantities, like magnetostriction and maximal tensile stress, show a smooth 
interpolation behavior between the two extreme cases. Additional reason for this hope is a very smooth interpolation 
behavior for magnetization and the penetration field that have already been calculated numerically for arbitrary 
thickness. For example, the normalized M{H) curves for all thickness/radius ratios differ by not more than 3% from 
some average M{H). It should be mentioned that our magnetostriction plots for long cylinder® and thin disk. Fig. 3, 
also show a remarkable similarity. 

In conclusion, we solve the magneto-elastic problem and find stress, strain and magnetostriction in a thin disk 
superconductor for the full cycle of applied field. The maximum tensile stress occurs approximately midway between 
the maximum field and the remanent state, see Fig. 5, and has quasi-parabolic dependence on the peak field, see 
Eq. (20). Qualitative behavior of the stress and strain components, and the normalized magnetostriction loop is 
similar to the case of long cylinder. 
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